Spectral density matrix of a single photon measured 
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We propose and demonstrate a method for measuring the spectral density matrix of a single 
photon pulse. The method is based on registering Hong-Ou-Mandel interference between photon 
to be measured and a pair of attenuated and suitably delayed laser pulses described by a known 
spectral amplitude. The density matrix is retrieved from a two-dimensional interferogram of coin- 
cidence counts. The method has been implemented for a type-I downconversion source, pumped by 
ultrashort laser pulses. The experimental results agree well with a theoretical model which takes 
into account the temporal as well as spatial effects in the source. 



Development of single photon sources brings the 
promise of implementing novel quantum-enhanced tech- 
nologies. In many applicationSjincluding quantum com- 
puting based on linear optics [l| photon sources are re- 
quired not only to deliver single light quanta, but also to 
supply them in a well defined mode. This is a necessary 
condition for quantum interference between independent 
sources [3| which is required in the above schemes. Also 
observation of three and more photon interference effects 
puts stringent requirements on the sources 01 • Besides, 
characterizing single photons is also interesting from the 
fundamental point of view. Historicahy photons were 
first described within the framework of quantum field the- 
ory (see [3| for a review) but more recently it was pointed 
out that a photon wave function can be introduced 
This kind of description, generalized by the introduction 
of the density matrix constructed out of projectors on 
the states with specific wavefunctions, seems to be the 
most elegant and effective theoretical tool for developing 
quantum-enhanced technologies 

Up till now measurements of polarization Q and spa- 
tial density matrix of a single photon were reported. 
The temporal characteristics of single photons were as- 
sessed only by verifying whether they interfere with other 
sources p, 1.1 (J| or between themselves In this Let- 

ter we propose and demonstrate a method for complete 
characterization of the temporal degree of freedom: a 
measurement of the spectral density matrix of a single 
photon. We show that the two-dimensional map of coin- 
cidence counts recorded as a function of delays between 
an unknown photon and a pair of weak reference pulses 
can be used to reconstruct the magnitude and the phase 
of the density matrix. We present a measurement for a 
type-I spontaneous down-conversion process in a bulk /3- 
barium borate (BBO) crystal, and compare the results 
of the reconstruction with theoretical predictions. 

When a single photon is launched into a singlemode 
fiber, the situation is significantly simplified since the 
spatial mode is well defined. If moreover polarization 
of the photons is fixed, the only remaining degree of free- 
dom is the spectral one. In this case a single photon 
component of the field can be described by the fohowing 



density operator: 

1 1 dujdijj'p{uj,Lo')a){uj)\()){Q\a{J) (1) 



where a(ijj) is an operator annihilating photon of fre- 
quency uj in the fiber, while p(uj, uj') is a density matrix 
of a single photon given in the spectral domain. 

Our method for measuring p{lu,uj') is based on the 
Hong-Ou-Mandel interference effect between the single 
photon to be characterized and a local oscihator (LO) 
pulse of known shape attenuated to a single photon level. 
The visibility of the two-photon interference dip is pro- 
portional to the overlap between the modes of interfering 
photons. Figuratively speaking, by modulating the local 
oscihator spectral amplitude and measuring the dip we 
can examine the single photon wavefunction from many 
directions. For a suitable class of LO pulses this suffices 
to retrieve the spectral density matrix of a single pho- 
ton. In a broad sense, our experiment is a single photon 
analog of the homodyne method for measuring quantum 
correlations within a hght pulse Indeed the den- 

sity matrix describes correlations within a single photon 
pulse. 

The method is presented in Fig. [TJ The first part is 
a Michelson interferometer which serves as a LO pulse 
modulator. The second part is the beamsplitter on which 




FIG. 1: (Color online) Scheme of the experiment. A weak 
laser pulses described by a spectral amplitude A{oj) enter a 
Michelon interferometer, comprising mirrors Ml and M2 and 
a 50/50 beamsplitter BSl, in which a double pulse of the local 
oscillator 0lo(i^) is produced. Then the LO interferes with 
an unknown photon described by the density operator p on 
a 50/50 beamsplitter BS2 and may give rise to a coincidence 
click of the detectors Dl and D2. 
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the Hong-Ou-Mandel interference occurs. When a mas- 
ter laser pulse described by a spectral amplitude function 
A{ui) enters the interferometer, it is split into two pulse- 
lets centered around delays ti and t2 determined by the 
length of the arms in the interferometer. If both pulselets 
overlap, they interfere on the beamsplitter and depend- 
ing on the precise value of their relative delay ^2 — ti the 
input pulse energy is directed alternatively towards the 
output of the interferometer or back to the laser. The 
normalized spectral amplitude of local oscillator pulses 
prepared this way reads: 



exp(— io^ii) -I- exp(— ia;i2) 



(2) 



where 5(^2 — ti) is a normaHzing function equal to the 
probability that a single photon will pass through the 
Michelson interferometer. The modulated LO pulses in- 
terfere with unknown single photons on a 50/50 beam- 
splitter BS2. If the master laser pulse contains a photon 
with probability I, while the source produces photon de- 
scribed by spectral density matrix p{lu,uj') with proba- 
bility / per pulse, we can calculate the probability of reg- 
istering coincidence with detectors of quantum efHciency 
ry to be 0|: 

Pc{tl,t2) - flSit2 - il)^ [1 - Q(tl, t2)] , (3) 

where we have assumed that each of the input 
states contains at most one photon and Q(ti,t2) — 
J J dwdu' (I>Iq{uj)p{ijJ,ijj')4>i,o{ijj') is an overlap between 
the local oscillator 0lo(^) and the density matrix 
p[ijj,Lo'). We can partially evaluate Q{ti,t2) using 
Eq. {21): 

, \ D{ti)+D{t2)+2Rep{hM) ... 
^^''^''^ 2S{t2-h) 

where Re stands for real part, while 
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is the inverse Fourier transform of a density matrix re- 
striced to the spectral domain of the master laser pulses 
and D{t) = p{t, t) is the function describing the standard 
Hong-Ou-Mandel interference dip. Inserting Eq. (0]) into 
Eq. {HI yields: 

Npc{h,t2) = 2S{t2 - h) - 2Rep{ti,t2) - D{ti) - D{t2) 

(6) 

where is a normalization factor combining the proba- 
bilities of detecting photon from the master laser I and 
the source /, as well as detector efficiency rj. At this 
point let us note that the density matrix p{u!,uj') is 
nonzero only near some point in the frequency space 



u! = Lo' = ujQ in case of narrowband photons, where 
loq is the central frequency of the single photons which 
we will assume to be equal to the central frequency of 
the master laser. Therefore Rep{ti,t2) oscillates like 
cos[a;o(t2 — ^i)]- Also S{t2 — ti) contains such oscilla- 
tions. On the other hand D{t) is a slowly varying func- 
tion of its argument. Therefore individual components 
of the right hand side of Eq. {U can be separated in the 
frequency domain. Let us apply the Fourier transform 
JJ dtidt2 exp(— ZLJifi -I- 2^2*2) ■ • ■ to both sides of Eq. {6]) 
and rearrange terms: 

A*{uJi)p{uJl,L02)A{uj2)+A{-UJi)p*{-UJi, -UJ2)A*{-UJ2) 
= (5(CJ1 - UJ2)\A{uJi + CJ2)P - Npc{ijJi,UJ2) + . . . (7) 

where by dots we have denoted Fourier transform of 
D{ti) and D{t2), while pc{uJi,uj2) is the Fourier trans- 
form of coincidence counts. Additionally we used the 
fact that is the Fourier transform of S{t). 

In the experiment we can directly measure Pc{ti, 12) as 
a function of ti and t2- In the case of narrowband sin- 
gle photons retrieving of the spectral density matrix runs 
the following way: first Fourier transform of Pc(ii,i2) is 
computed. A region in the frequency space where con- 
tribution from S(t2 — ti) and p{t\,t2) lies is separated. 
Next, a large contribution from S{t2 — h) is calculated 
and subsequently subtracted by measuring Pc(ti,t2) for 
large ti and t2 where p(ti,t2) is zero. This way we 
obtain the Fourier transform of p{ti,t2), which equals 
A* {uji)p{ijJi,uj2)A{ijj2). It is divided by spectral ampli- 
tude of master laser pulses A{lo), which is measured sep- 
arately and finally p{uji,uj2) is found. Note that the last 
step is well defined only for single photons of bandwidth 
narrower than that of the master laser. 

To illustrate the above dry formulas a typical coin- 
cidence pattern and its Fourier transform are depicted 
in Fig. [21 The diagonal fringes in Fig. [2ta) come from 
S{t2 — ti), the vertical and horizontal stripes from D{ti) 
and Dlt2) while the most interesting term p(ti,t2) con- 
tributes only in the very center of the picture. It is easier 
identified in Fourier transform plot in Fig. [2jb) where 
a diagonal cloud corresponds to ^*(cji)p(ciJi, [^2)^(^2), 
D{ti) and D{t2) contribute a cross arround zero fre- 
quency while Fourier transform of S{t2 — ti) is seen as a 
ridge along oji = uj2- 

In the above derivations we have assumed interference 
between single photons. In our experiment, in fact we 
interfere a weak coherent state with a multimode ther- 
mal state. Whereas the visibility of such interference can 
be exactly calculated using the semiclassical theory 
two photon terms of the coherent state and the thermal 
state contribute only towards a constant background of 
coincidence counts. However, the shape of the interfer- 
ence pattern does characterize the single photon compo- 
nent of the signal field. Its density matrix can still be 
retrieved in the way described above. 
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FIG. 2: (a) A typical coincidence interferogram as a func- 
tion of delays ti and t2 and (b) its Fourier transform. The 
black parallelogram in (a) defines the scan range, which with 
suitable sampling density yields the region of interest in the 
frequency domain, outlined in (b). 



Our experimental setup is depicted in the Fig. [31 The 
master laser (RegA 9000 from Coherent) produces a train 
of 165 fs FWHM long pulses at a 300kHz repetition rate 
centered at 774 nm, of which we use 300 mW average 
power. Most of the energy goes to the second harmonic 
generator, based on a 1 mm thick BBO crystal cut for 



BS FL XSH IL 




FIG. 3: (Color online) The experimental setup. BS, beam- 
splitter; FL, focusing lens; XSH, BBO crystal for generation of 
the second harmonic; IL, imaging lens; DM dichroic mirrors; 
BG, blue glass filter; X, downconversion crystal; IF, interfer- 
ence filter; FC, fiber coupling stage, HWP, half waveplate; P, 
polarizer; ND, neutral density filter; FPC, fiber polarization 
controler; Dl and D2, single photon counting modules. 



0.8 



0.6 



0.4 



0.2 







/ 


/ 1 
/ / 

/ 

1 


N \ 
\ \ 
\ V 


/ 

/' . 


1 f 
Iff 
ll if 

— V 1. 







760 765 770 775 780 785 
X [nm] 



-12 



FIG. 4: Spectral intensity |A(aj)|^ (solid line) and phase 
arg A(tj) (dashed line) of the master laser pulse retrieved us- 
ing FROG. 



type-I process. Ultraviolet pulses produced this way have 
1.3nm bandwidth and 30 mW average power. They are 
filtered out of fundamental using a pair of dichroic mir- 
rors DM and a color glass filter BG (Shott BG39), and 
imaged using 20cm focal length lens IL on a downcover- 
sion crystal X, where they form a spot measured to be 
155 fim in diameter. The crystal X is a 1 mm thick BBO 
crystal cut a 29.7° to the optic axis, and oriented for max- 
imum source intensity. A portion of down-converted light 
propagating at an angle of 2.5° to the pump beam passes 
through a lOnm interference filter centered at 774. 5nm 
and is coupled into a single mode fiber. This defines the 
spatial mode in which the down-conversion is observed 
[131 . About 4% of energy of master laser pulses is re- 
fiected towards an LO preparation arm. The pulses first 
go through a half wave plate HWP and a polarizer P 
allowing for fine control of the energy. Then they are de- 
layed in a computer-controlled delay fine. Next the pulses 
enter a Michelson interferometer in which one mirror 
moves on a precision computer controlled stage allowing 
for generation of double pulses with a well-defined tem- 
poral separation At = t2 — ti. Finally the double pulses 
are attenuated to contain less than 0.1 photon on average 
and coupled to a single mode fiber, where their polariza- 
tion is adjusted using a fiber polarization controller FPC 
to match the polarization of the photon coming from the 
downconversion source. Both the downconversion and 
local oscillator photons interfere in a 50/50 singlemode 
fiber coupler and are detected using single photon count- 
ing modules SPCM (PerkinElmer SPCM-AQR-14-FC) 
connected to fast coincidence counting electronics (suit- 
ably programmed Virtex4 protype board ML403 from 
XiHnx) detecting events in coincidence with master laser 
pulses. 

For calculating the actual density matrix of unknown 
single photons, a characterization of master laser pulses 
was necessary. This was accomplished using Frequency- 
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FIG. 5: (Color online) Contour plot of the measured spec- 
tral density matrix p{u>,uj') as a function of wavelengths 
A = 2-kc/lij (solid lines) and a theoretical prediction for this 
quantity (dashed red lines). The contour were drawn at 0.75, 
0.5 and 0.25 of the maximum values, the outermost encircles 
7x5 experimental data points. 



Resolved Optical Gating (FROG) technique fiH^. Re- 
trieved spectral intensity |A(a;)p and phase aig A{uj) are 
plotted in Fig. |4l Precise knovirledge of the master laser 
pulses allowed us also to calculate the second harmonic 
pulse using perturbative approach. 

The complete measurement consisted in a series of 6 
scans of a rectangular grid depicted in Fig. [2]Ja) spanned 
by 4000x25 points, where the latter number refers to the 
direction along the fringes. The corresponding mesh was 
0.233 fs X 66 fs, and coincidences were counted for 80 ms 
at each point. The reconstructed spectral density ma- 
trix of a single photon is depicted in Fig. [H We compare 
it with theoretical calculations plotted with dashed lines 
in the same figure. The theoretical model used in these 
calculations assumed the exact phase matching function 
of the nonlinear crystal and the ultraviolet pump pulse 
shape computed from the measured A{uj). The trans- 
verse components of the wave vectors for the pump and 
down-converted beams were treated in the paraxial ap- 
proximation. The spectral density matrix was calculated 
for coherent superpositions of plane- wave components of 
the down-conversion light that add up to localized spatial 
modes defined by the collecting optics and single-mode 
fiber. The other photon from the source, which remains 
undetected, was traced out assuming that it can prop- 
agate at any direction and have any frequency that is 
consistent with the conservation of energy and perpen- 
dicular momentum in the downconversion crystal. As 
seen in Fig. [5] the theoretical calculations predict more 



pronounced correlations i.e. smaller width along the an- 
tidiagonal p{uJo + S,ujo — S) than was actually measured. 
We attribute this discrepancy to a difference between the 
actual ultraviolet pump pulse shape and the one calcu- 
lated from A{u!). Also the tips of the density matrix are 
measured with reduced accuracy, since in that regions the 
raw experimental result is divided by a relatively small 
master laser spectral intensity A{uj) which amplifies er- 
rors. The theoretical model predicts the phase of the den- 
sity matrix to be smaller than tt/SO in the region bounded 
by the contour at 0.25 maximum. The measured phase 
in this region is smaller than tt/IO and varies randomly 
from point to point. 

In summary, we proposed and demonstrated a method 
for measuring the spectral density matrix of a single 
photon component of the electromagnetic field in a sin- 
glemode fiber. The method is based on two photon 
interference and is thus limited to the spectral range 
where known reference pulses are available, however it 
allows for retrieving both amplitude as well as phase of 
the density matrix. We have applied this method to a 
downconversion-based source of single photons and found 
that measured density matrix agrees with theoretical pre- 
dictions. 
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